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Abstract 



Two planar supersymmetric quantum mechanical systems built around the quantum 
integrable Kepler/Coulomb and Euler/Coulomb problems are analyzed in depth. The 
supersymmetric spectra of both systems are unveiled, profiting from symmetry operators 
not related to invariance with respect to rotations. It is shown analytically how the first 
Q^ problem arises at the limit of zero distance between the centers of the second problem. It 

appears that the supersymmetric modified Euler/Coulomb problem is a quasi-isospectral 
deformation of the supersymmetric Kepler/ Coulomb problem. 
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1 Introduction 

During the last forty years a very interesting jump from symmetry to supersymmetry has taken 
place, determining theoretic particle spectra in quantum field theories with extremely appealing 
characteristics, see e.g. [I]. Unlike many quantum field theoretical models, the supersymmetric 
systems are frequently amenable to non-perturbative treatments, see e.g. [2], but the main fea- 
ture is that fermions and bosons are jointly assembled in multiplets, a fact, although suggestive, 
that has not yet experimentally confirmed. Thus, mechanisms of spontaneous supersymmetry 
breaking must be investigated in the search for explanations of the apparent lack of super- 
symmetry in nature. In a series of papers, Witten, [3], |1], and [5], proposed the analysis 
of this phenomenon in the simplest possible setting: supersymmetric quantum mechanics. A 
new area of research in quantum mechanics was born, with far-reaching consequences both in 
mathematics and physics. The relation between the Dirac operator in electromagnetic and/or 
gravitational fields - the supercharge - with the Klein-Gordon operator - the supersymmetric 
Hamiltonian - provided a guide for the building of supersymmetric quantum mechanical sys- 
tems. The factorization method of identifying the spectra of Schrodinger operators by means 
of first-order differential operators, see [6] for a review, is another antecedent of supersymmet- 
ric quantum mechanics that can also be traced back to the 19th century through the Darboux 
transform [SO]- In its modern version, supersymmetric quantum mechanics prompted the study 
of many one-dimensional systems from a physical point of view. A good deal of this work can 
be found in References [7], [8], [9], [in]. Several examples of this structure with emphasis in 
the semi-classical behavior of non-harmonic oscillators have been worked out in [Tl] . 

The formalism of physical supersymmetric systems with more than one bosonic/fermionic 
pairs of degrees of freedom was first developed by Andrianov, loffe and coworkers in a se- 
ries of papers, [12], [13], pubhshed in the eighties. Factorability, even though essential in N- 
dimensional SUSY quantum mechanics, is not so effective as compared with the one-dimensional 
situation. Some degree of separability is also necessary to achieve analytical results. For this 
reason we started a research program in the two-dimensional supersymmetric classical me- 
chanics of Liouville systems, [16]; i.e., those systems separable in elliptic, polar, parabolic, or 
Cartesian coordinates, see papers [TTj and [IB]. We therefore follow this path in the quantum 
domain for Type I Liouville models in [TH] . 

Nevertheless, the authors from Saint Petersburg University mentioned above considered 
from the earlier eighties higher-than-one-dimensional SUSY quantum mechanics from the point 
of view of the factorization of N-dimensional quantum systems, [H], [15]. loffe et al. also stud- 



ied the interplay between supersymmetry and integrability in quantum and classical settings in 
other types of model in References [20], (2T], [22]. In these papers, a new structure was intro- 
duced |2Hj: second-order (and higher-order) supercharges provided intertwined scalar Hamil- 
tonians even in the two-dimensional (and higher-dimensional) case, see the review papers [23] 
and [21]. This higher-order SUSY algebra allows for new forms of non-conventional separability 
in two dimensions. There are two possibilities: (1) a similarity transformation performs the 
separation of variables in the supercharges and some eigen-functions (partial solvability) can be 
found, see [23], [2S]- (2) One of the two intertwined Hamiltonian allows for exact separability: 
the spectrum of the other is known, [27], [25] . 

Our purpose in this paper is to describe planar supersymmetric systems - two bosonic/fer- 
mionic pairs of degrees of freedom - such that the Bose-Bose and Fermi-Fermi (scalar) Hamil- 
tonians will be separable. In Reference [31] Eisenhart classified all the quantum systems with 
separable Schodinger equations in Cartesian, polar, parabolic, and elliptic coordinates. We 
shall address two planar supersymmetric separable systems, one in polar, the other in elliptic 
coordinates. 

We shall first consider the planar supersymmetric Kepler/Coulomb problem showing the 
separability in polar coordinates. We strongly rely on the work by Wipf et al in papers [12], 
|43j where they solve this problem in any dimension D by finding a supersymmetric matrix 
Runge-Lenz vector and describing algebraically the spectral problem in terms of the irreducible 
representations of SO{D + 1). Instead we shall attack the spectral problem in the Bose-Bose 
sector, finding the bound state energies in terms of the Casimir eigenvalues of the irreducible 
representations of SO(3), whereas the eigenfunctions are generalized Laguerre polynomials. The 
scattering eigenfunctions in this sector, as well as in the Fermi-Fermi sector, are generic conflu- 
ent hypergeometric functions given in terms of infinite series. By acting with the supercharges, 
we provide in turn both the bound state and scattering eigenfunctions in the Fermi-Bose Sector. 
We remark that, following a previous work on the supersymmetric classical Kepler/Coulomb 
problem, [H], Heumann chose another superpotential [10] leading to a supersymmetric quan- 
tum mechanical system where the Runge-Lenz vector is no longer an invariant even in the 
Bose-Bose and Fermi-Fermi sectors. 

In the second half of the paper we study a supersymmetric quantum mechanical system 
built from the classical Euler problem: a light particle moving in the gravitational field created 
by two fixed Newtonian centers of force restricted to the plane of the centers [15]. Besides 
Euler, this system attracted the interest of investigators of stature such as Lagrange, Jacobi, 
Liouville, Darboux and others, see [IS] to read a brief history of the subject, on a double front: 
1) because of the potential applications in celestial mechanics, e.g., as an intermediate step 
in the three-body problem. 2) Because the Euler problem was a playground where the ideas 
of integrability, curvilinear coordinates, Hamilton- Jacobi separability, of such importance in 
classical dynamics, were tested. All this was imported by Pauli [IH] to the quantum domain 
in his research on the spectrum of the H2 hydrogen ion molecule. A Chapter of Pauling and 
Wilson's book on Quantum mechanics [17] is devoted to the developments in this quantum 
problem up to the mid thirties of the past century. 

We shall address a supersymmetric quantum mechanical system such that the scalar Hamil- 
tonians in the Bose-Bose and Fermi-Fermi sectors are related to the quantum mechanical Eu- 
ler/Pauli Hamiltonian. We are guided by the separability of the Schrodinger equation in elliptic 
coordinates: half the sum and half the difference of the distance to the centers. This is the 
main property of the Euler-Pauli Hamiltonian allowing for its integrability. We choose our 
scalar Hamiltonians fulfilling this property but supersymmetry requires the energy to be non- 
negative. We are forced to add a "classical" piece to the EP potential energy that pushes the 
ground state energy to zero. All this is achieved by the choice of a superpotential inspired 
in the loffe/Wipf et al superpotential for the supersymmetric Kepler/Coulomb problem (our 



superpotential tends to the ABI/KLPW superpotential when the two centers collapse). In 
Reference [52], however, we explored other possibilities in comparison with this superpotential. 

A double change of variables to one-half of arccosh and arccos of the elliptic variables trans- 
forms the separated spectral problem into systems of entangled Razavy |1H] and Whittaker-Hill 
|49j (three-term Hill, Razavy trigonometric) equations. The two equations in each system are 
in principle independent, but they are entangled because their parameters are determined from 
the integration constants -the energy and the separation constant- that allow one to formulate 
the Hamiltonian spectral problem in terms of two separated ODE's. To the best of our knowl- 
edge, this change of variables was formulated for the first time in [32] . We shall profit from the 
fact that for certain values of the parameters either the Razavy or the Whittaker-Hill equation 
are algebraic quasi-exactly solvable systems q More specifically, the bound states of our sys- 
tem arise from values of the energy and the separation constant that lead to Razavy equations 
belonging to this class of algebraic QES potentials, see [38] . The algebraic QES potentials have 
corresponding quantum Hamiltonians that are elements of the enveloping algebra of a finite 
dimensional Lie algebra (Lie §L(2,]R) in many cases), admitting an invariant finite module of 
smooth functions as irreducible representations [37j. These potentials were first studied by 
Turbiner [33] and then completely classified by Gonzalez-Lopez, Kamran, and Olver [M]-|35]. 
The associated (weakly) orthogonal polynomials were analyzed in full generality in [36] and all 
this machinery was applied to study the Razavy trigonometric potential in [39] . 

Unlike in the non supersymmetric case, see [5l]-[H5], the Whittaker-Hill equations unfortu- 
nately are not QES for the values of the parameters for which the Razavy equations are QES 
in the supersymmetric case. Thus, we can only give the solutions in the form of infinite series 
following the theory of Hill equations, see e.g. [56]- [57]. The case of two centers of the same 
strength is exceptional: instead of Whittaker-Hill equations we encounter Mathieu equations. 
The solutions can be given analytically in terms of Mathieu Cosine and Sine functions F|53j 
and it can be explicitly checked that at the limit where the two centers collapse the one-center 
wave functions are recovered. 

The organization of the paper is as follows: after this long Introduction in the second 
Section §.2 we settle down to the framework of A/" = 2 supersymmetric Quantum Mechanics 
for systems of two degrees of freedom. Because each degree of freedom can be labeled either as 
Bosonic or Fermionic, we have 2^^ = 4 types of state. The Clifford algebra of M"^ describes this 
situation perfectly and helps us to define the supercharges, the Hamiltonian structure, and the 
Hilbert space of states. Section §.3 is fully devoted to discussing the planar supersymmetric 
Kepler/ Coulomb problem. The Bose-Bose bound state eigenfunctions arise as irreducible repre- 
sentations of the dynamical §0(3) symmetry associated with the Runge-Lenz vector whereas 
Bose-Fermi strictly positive bound states are obtained via the action of the Q^ operator on 
the BB bound states (except the zero mode). The scattering states are also describe to unveil 
the whole spectrum of the supersymmetric Kepler/ Coulomb problem. A direct analysis of 
the spectrum of the 4x4 matrix Hamiltonian and symmetry operators is also included. In 
Section §.4 we address the same program in the supersymmetric modified Euler/Pauli two- 
center system. The bound state wave functions come from polynomial x exponential solutions 
of (infinite) Razavy equations multiplied by power series solutions of related Whittaker-Hill 
equations. In the case of two centers of the same strength, the WH equations are replaced by 
Mathieu equations. In Section §.5 we show how the Kepler/Coulomb spectrum reappears at 
the d = limit of the two centers. Finally, we offer a final Section with further comments on 
interesting generalizations of these classicalljjj integrable models. 



^In the non supersymmetric problem Demkov in |54j has shown that there are exceptional values of the 
energy and the separation constant for which both equations in the pair are QES providing finite solutions. 
^All the conventions on special functions throughout the paper will follow Reference [53] , 
■^Here we use the word classical in a non physical sense, i.e., classical does not refer to a class (versus quantal) 



2 J\f = 2 supersymmetric planar quantum systems 

2.1 J\f = 2 two-dimensional SUSY quantum mechanics 

In this type of quantum mechanical systems there are two pairs of canonically conjugated 
Bosonic operators - giving the position and momentum of the particle - that we choose in 
coordinate representation: 

d 
Pk = -ih- — , Xk = Xk, [xk,Pj] = ihSkj . 

There are also two pairs of Fermionic operators -of physical dimensions: [ipk] = M~2 - taking 
care of the Fermionic degrees of freedom of the system. The Fermi operators satisfy anti- 
commutation relations of the form: 

{4, ^J = = {^l ijj} , {4, i^h = l^hi , fc, / = 1, 2 , (1) 

showing that one operator is canonically conjugated to its adjoint operator. 

The Fermionic Fock space is built from the vacuum state: ipk\0) = ,VA; = 1,2, the two 
degrees of freedom being in Bosonic states because |0) is an eigenstate of zero eigenvalue of the 
Fermi number operator N = X]fc=i Yk^k'- N\0) = 0|0). The creation operators acting on |0) 
bring the system into one-particle states ■?/'|,|0) = |lfc), where one of the two degrees of freedom 
becomes Fermionic: N\lk) = |lfc)- The two-particle state - the two degrees of freedom in 
Fermionic states A^|lil2) = 2|lil2)- are then obtained in a dual way related by Fermi statistics: 

V'||ll) = |lll2) = -^I|l2) = -|l2ll). 

The ortho-normality relations 

(0|0) = (li|li) = (l2|l2) = (lll2|lll2)=l 

(li|0) = (I2IO) = (I2II1) = (I1I2IO) = (I2I1II1) = (I1I2II2) = 

allow us to write the more general state in this finite Fermionic Fock space J-" = Tq (B J^i (B J^2 
in the form: 

2 

l/) = /o|0) + $^/lfc|U)+/2|lll2) , fo,fik,f2eC . 

k=l 

The supersymmetric space of states is the direct product of J-" with the Hilbert space L^(]R^): 
STi = 1-1® J^ = L'^iM?) (g) C^ = SHq © SHi © S1-L2- The supersymmetric wave functions read: 



|^(xi,X2)) = /o(a;i,a;2)|0) + ^/ifc(xi,X2)|lfc) + /2(a;i,a;2)|lil2) 



fc=i 



/o(a;i,X2),/ifc(xi,X2),/2(a;i,X2) G L 



2/Tn)2\ 



The standard procedure for introducing supersymmetric dynamics in this setup runs as 
follows: One first defines the supercharge^^ 

V OXk OXkJ V 9Xk OXkJ 



of physical phenomena. 

^By [O] we shall denote the physical dimensions of the observable O. 



where 



W{xuX2) : 



— )■ 



[W] = ML^T-^ 



is a real function called the superpotential. The supercharges are thus nilpotent first-order 
differential operators Q^^ = Q = Q^ that move states between the different Fermionic sectors of 
the space of states: 

O riQ ^ C> ri\ ^ o rL2 

Q Q 



The super-Hamiltonian is defined to be 

1 



H 



1 



^{Q,Q^] = ^[QQ^ + Q^Q 



which implies: [if, Q^] = [H, Q] = 0. Therefore, the transformations generated by Q and Q^ 
are symmetries - called supersymmetries - of the dynamical system and the supercharges are 
themselves constants of motion. Because [H, N] = 0, there is a Z2-grading of the dynamics 
given by the Klein operator F = (—1)^: 



F|0) = |0) , F|lil2) = |lil2) , F|li) = -|li) , F|l2) 



II2) 



In the classification of Supersymmetric Quantum Mechanics given by Kibler et al. in [50] our 
formalism ranks in the class of a complex super-charge, Q, with an involution operator F"^ = I. 
It is also shown in Reference [50] that it is equivalent to another supersymmetric system with 
two real supercharges: this is the reason for the TV = 2 in the title. 

2.2 Clifford algebra representation 

In order to skip abstract ket/bra Dirac algebra we represent the Fermi operators by means of 
the generators of the Clifford algebra of M^: 



V'l 
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2^/m 
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/ 1 \ 
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0-1 



\ / 



One can check that this is a minimal realization of the Fermionic anticommutation rules M 
and the Fermionic Fock space becomes the space of four-component Euclidean spinors with 
basis: 
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The supercharges are 4 x 4-matrices of differential operators 
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f Di D2 \ 
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\ -D2 i^i / 



dW 
da 

differential operators 



where Dk = hdk + |^ and Dk 



hdk — 1^. The super-Hamiltonian is also a 4 x 4-niatrix of 



H = Ha 



fc=l 1=1 



dxkdxi 



i^m 



f Ho \ 

Hl^ Hl^ 

Hf^ Hf 

\ /J2 / 



with a block-diagonal structure inherited from the eigen-spaces of the Fermi number operator: 



N = V-lV-i + V-a^ 



Thus, in the F = +1 eigen-sectors of STi the Hamiltonian act by means of the scalar ordinary 
Schrodinger operators: 
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S-H2 2m 
In Sl-Li, however, the super-Hamiltonian reduces to the 2 x 2-matrix Schrodinger operator: 



Hi = H 



SHi 



Hi' Hl^ 
Hf Hf 



H. 



fi. Q2 



diW - 
-d2diW Ho 



-did2W 



S Q2 
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We see that all the interactions come from the gradient and the second-order partial derivatives 
of the superpotential. 

3 The planar quantum Kepler/Coulomb problem and 
supersymmetry 

Our first goal in this survey is the development of this formalism encompassing the Hamiltonian 
of the Kepler-Coulomb problem. 



3.1 The quantum Kepler/ Coulomb Hamiltonian 

We recall that the Kepler/Coulomb Hamiltonian describing the quantum dynamics of one- 
electron atoms is: 



K = — ^A-- ,a>0 , \a] = ML'T 
2m r 



3rp-2 



We re-scale positions and momenta to new variables x^ — ?■ -^x^, Pk -^ mapk with dimensions 
of [xk\ = M'^L''T~'^ and [pk] = M~'L~'^T. By this token we see that the parameters m (particle 
mass) and a^ (strength of the coupling) factor out in the new Hamiltonian 



K -)■ ma^K = ma'^ ( A 



1 

- A -- 
2 r 

and their only physical role is to set the energy scale 



[K] = M-^L-^T 



It is well known that this problem is superintegrable : The angular momentum -one scalar 
in the plane- 

L = —in Xi— X2- 



8x2 dxi 

and the Runge-Lenz vector -two components in the plane- 

^1 = 2 [P^L + Lp2J -— , ^2 = "2 [PiL + Lpij - — 

[L, ii] = ifiA2 , [L, is] = -ifiK 

both commute with K: 

[i^,L] = [K,ii] = [K,i2]=0 . 

We remark that in our variables the physical dimensions of these operators are [L] = ML'^T~^, 
[Ai] = [A2] = 1 and recall that they close the SO(3) Lie algebra in the space of negative energy 
(bound states) eigen-functions of KipE = EipE, E < 0: 

[A,, A2] = -m fd±A -l)^ = -'^^^KL , Ml = ^^A, , M2 = yi^^ 

M3 = L , [Ma, Mfe] = theabcM, . 

Moreover, because the SO (3) Casimir operator is 

C^ = Ml + Ml + Ml = -X {Al + Al 

and we have 

Al + Al = 2k{L^ + ^\+l 

the Hamiltonian is given in terms of C^ 



2 (72 + ^ 

such that the SO (3) symmetry is not Noetherian but a dynamical symmetry. One finds imme- 
diately the bound state eigenvalues 

^' = 4 ■ (2jTT)5 • '^l^"" ' 

which must be multiplied by ma'^ to find the physical bound state energies. The bound state 
eigenfunctions, the SO (3) irreducible representations, will be given in the next subsection. 

3.2 The supersymmetric Kepler/ Coulomb Hamiltonian 

The superpotential proposed by loffe et al in [15] and Wipf et al (independently) in |12] to 
build the supersymmetric version of the Kepler/Coulomb problem isj 

W(xi X2) = --r = -- fil^^^ — = --^ ^'^ ^ 2 / XfcXA 

h h^ ' dxk h r ' dxkdxi hr \ r^ J 



5 



In [30] Heumann proposed another superpotential which spoils the Runge-Lenz vector conservation. 



We also re-scale the Fermi operators ipk — ;• -^"^fe 
supercharge: 

d 2 



where 
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[z/'jt] = 1 to define the Kepler-Coulomb 
[Q] = M-^L-'^T 
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"tpk 



9 







is', 9} 
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2 dxk r 

is a "hedgehog" projection of the spin variables over the M^-plane. Explicitly. 
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The supersymmetric Kepler/ Coulomb Hamiltonian reads: 



H 



h^ 



A 



2 
1)? 



L X 



X=\gr9^]=lA-2N + rg^9 



2 n^y r 

whereas the scalar Schrodinger operators in the Bosonic sectors are: 



h^ 2 



K + 



h? 



2 n^ r 



r> 2 2 



(2) 



Thus, ^0 is exactly the Kepler/ Coulomb Hamiltonian plus a constant needed to set to zero the 
energy of the Bosonic ground state (zero mode); recall that supersymmetry forbids negative 
energy eigen-states. if 2, however, is also (modulo a constant) the Kepler/ Coulomb Hamiltonian 
for a particle of opposite electric charge, say a positron. The force is repulsive and there will 
only be scattering states. 

The matrix Schrodinger operator - already given in [15] circa 1984 - acting in the two- 
dimensional sub-space of the Fermionic Fock space such that A^ = 1 is: 



Hi 



2x1X2 



1 ^2 



2x1X2 

2 ^ -h ^2 -h ^3 



3.3 A^ = bound state eigenfunctions of H 

The A^ = bound state eigenfunctions of H are exactly the eigenfunctions of ifo, which are 
the same as the bound state eigenfunctions of K with displaced eigenvalues: 



H,^) 
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n-^ 2 h\j{j + i) 
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Vn2 2 h\j + \y 



^) 
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n 1 
J = 7T e - 



N , ^o^f = ma'Ef^pf 



2ma 



1- 



(2j + l)V ' 



(0) 



Because of the dynamical SO (3) symmetry of Hq, the bound state eigen-functions, which are 
degenerated in energy, form irreducible representations characterized by two integer or half- 
integer numbers, j and m, providing the eigenvalues of the Casimir operator and M3 in common 
eigen-kets: 



C'^\i]m) = h^j{j + l)\j]m) , Mslj-^m) =hm\j]m) , m : -j,-j + 1,- ■ ■ ,j - l,j 



Using polar coordinates r = +A/xf + x|, if = arctan— in coordinate representation the eigen- 
wave functions are of the form: ^jJi{r, (p) = (r; (p\j; m). It is not difficuh to identify the highest 

weight eigen-wave functions in each irreducible representation. Let M+ = Mi + iM2 be the 
up-stairs ladder operator 

M+ = hi J + -] A+ = K" [ J + -] e'^P li- 



1) \ 2/ [ drdif r dif'^ 2 dr 2r dcp K^ 

which annihilates the highest weight state: 

MH_V^gHr,y.)=0 , V^g)(r,y.) = (r;¥.|j;j) = /,(r)e^^''^ . 
This first-order ODE is easily integrated 

and the normalized highest weight wave functions are: 

^S?(r,y,)=2./?- , ^' -e-^-^-e-t , « = —ilL- . (3) 

The down-stairs ladder operator M_ = Mi — iM2 in polar coordinates reads: 

\ 2/ \ 9r(9v9 r (9^9^ 2 9r 2r d(p h? 

From the Lie algebra we see that 

il)f^ oc M^j'^ipfj , m = j,m = j - 1,- ■ ■ ,m = -j + l,m = -j , M^'^^V'jj = 

and from the ansatz ilj-J^{r,(f)) = N^~^^^r^"^^Pj_\jji\{r)e 1^(2^+1)6*'"'^, where the iV-^"!™-! are nor- 
malization constants, the recurrence relations 



Pj-\m\{r) 



PH.i„H.,-?(^±M^; 



P,-»-M . £t4±illp' 



j-|m|-l 



T 



follow. Therefore, the N = bound state eigen-functions are: 



In Q we have the Kummer confluent hypergeometric functions iFi for those values of a and b 
such that the series 

^r , 1 ^ ia)k z^ ,x r(a + l) „, r(6+l) 

truncate to generalized Laguerre polynomials: 

In sum, the A^ = bound state eigen-functions of the planar supersymmetric Kepler-Coulomb 
problem are organized as (degenerated in energy multiplets) irreducible representations of Sd 
in L^(M"'' X §^) rather than in L^(§^) (spherical harmonics). 

10 



3.3.1 Ortho-normality and lower energy levels 

It is easy to clieck that the following ortho-normality relations hold: 
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Note that in the case of one integer ji and one half-integer j2 pairing it is necessary to integrate 
the p variable over Air because of the double-valued representation. 

We now offer two Tables with the lower energy multiplets, their 2D (cross-sections) and 3D 
plots: 



Energy 



Eigen-function 



4") = 



,,(0) 



noiT^^)^w\i e « 



M°'= ' 



2h? 






E. 



(0) _ j^ 



V'i°^(r,^) 



8 



m'^y^ 



r e 3"^ e 



'<p[°-iir,f) = 



<\r,^)-i^Ji. i^f^'-^r) e-^ 



9h'^V37T 



r e 31-' e 



E' 



(0) _ j^ 
" 8R2 



r3/2 



f^/^ g 2h'^ e 2^ 






If 

V'f L = 5^ ri/2 {r - 2fi2) e-^ e'^^ 

„3/2 p_~^^ „— 2|t^ 



2 2 

^^3 



4ft^y37r 



p(0) 
^2 



771^"^ _ 48 



V'2?(^^) 



16 



-^22 viiry 125 fts 7l5^ 



r^ e~ss^ e^^'^ 



^20^ (^' ^) = Tsfl^ \/^ (25 ft* - 40 fi2 r + 8 r2) 






■02-2 (^>'/') 



16 



125rt«Vl57r 



r e 5fi^ e 



-2i^ 



3.4 A^ = 1 bound state eigenfunctions of H 

The supersymmetric partner eigen-states belonging to STii with the same energies are of the 
form: 



E) 



(1) 



h^ 



(2jTT?)"^°'|)'^'>°' 



/ 



■,(1)^^ ^A - ftnt./,(o), 



^P)rniXl,X2) = hQ^llj)J^{Xi, X2) 



V 









\ 
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Table 1: 3D Plots and cross-sections of the probability densities: A^ = 0. 



Probability 
density 



.,(0) 



moixi,X2)\' 



2-^2 



',(0) 



(0), 



|V'i±i(a;i,X2)r \iIjIq\xi,X2)\ 



h = i 





tW 





h= 1 




Therefore, 

■,(1) 



because Q'?/'j^(r, 99) = and H ipjj^{r^(p) = ^j ^ ipjj^{r^ip). Integration over M^ gives 
I Jo Jo Jo 

_ o fc2 771(1) r r 



i2 "ilJ2"mim2 



The normalized eigenspinors V'jM'"' v) ^^ / "^ ,^. "ipjmi''"^ ^) 



e~5 



— 7-1- 777, 

COSV9 < — — . iFi [-j + |m| + l,2|m| + 2,d + j iFi [-j + |m|, 2|m| + 1, 

[ 2\m\ + 1 

in 

H e""^ iFi [-j + \m\,2\m\ + 1, 

u 



u 



,u\ 



i^''lm'^\r,V>) =i 



Kj + HV- 



u 



\m\ 



smip 



m 



n? \j (2j + l)3((2j + 1)2 - l)(j - |m|)!7r (2|m|)! 
—j + \m 



e*'"^e-2 



2|m| + 1 



iFi [— j + \m\ + 1, 2|m| + 2,u] + j iFi [—j + \m\, 2\m\ + l,u] 



+i — e"**^ iFi [-j + \m\, 2\m\ + 1, u] 



u 



satisfy the spectral condition H ijjjj^{r,(p) = E- ijjjj^{r,(p) = Hi ^jjjj^{r,(p) and form an or- 
thonormal basis in SHi. 
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Specifically, these two-component wave functions are linear combinations of two contiguous 
generalized Laguerre polynomials. The reason is that Q does not commute with the generators 
of the SO(3) symmetry: [M^, Q^] 7^ 0, Va = 1, 2, 3. Therefore, the Q^ action does not respect the 
§0(3) irreducible representations. Nevertheless, the spinorial wave functions are characterized 
by the quantum numbers j and m, although the degenerated multiplets do not form irreducible 
representations of SO(3). We show next the lower spinorial probability densities: 



Table 2: 3D Plots and cross-sections of the probability densities: A^ = 1. 



Probability 
density 



|V'S(a;i,X2)P |V'2i(a;i,X2)P r|V'l^ji(xi,X2)P iV'l + sla^bXz)! 



±1 
2-^2 



2-^2 



h=i 







h=i 






3.5 Scattering states and supersymmetric Hodge spectral decompo- 
sition 



On positive energy eigen-functions of the Kepler-Coulomb Hamiltonian K the normalized com- 
ponents of the Runge-Lenz vector and the angular momentum close the §0(2, 1) Lie algebra: 



Ml 



:Al , M2 



--Ao 



M, 



V2E V2E 

[Ml, M2] = -ihMs , [Ms, Ml] = ^hM2 , [M2, M3] = ^^Mi 
To search for the scattering wave functions in STio, the eigenfunctions of Hq 



2 



2 



(0) 



<i,(xi,X2)=i?(°VS)(xi,X2) , i?(°)> 



h^ 



we profit from the fact that the spectral problem is separable into polar coordinates. The ansatz 
'i/'Lo)(r, (/)) = /Lo)('")e*'"'^ leads to the ordinary differential equation 



2 ^(0) 
E(0) 



u J 



dr"^ 



+ 



Idf^ 



(0) 

£;(o) 



r dr 



^2 Jew y) + 



h? \r 



1 



+ ^(0) 



h^ 



f 



(0) I 



Now defining the non-dimensional variable u 
solutions 



2^2\^-Em\ 



r we find the Bosonic scattering 



/ 



(0) , 



N{E^''^)e-'^u^"'KFi 



\7n\ 



1 
2 



fi^/2\^ 



E(o)| 



:, 1 + 2\m\,iu 



(5) 
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in terms of Kuinmer confluent hypergeometric functions. 

Simili modo, we identify the scattering wave functions in S'H2-, the eigenf unctions of H2'- 



^+4 + -VFi)(^i'^2)=i?('Vl%(a:i,X2) , E^^^ 






dr^ ^ r dr r^ ^^"'^^ ^ h^ \ r ^ h^ ^ ^'''^ ' ~ 



f%{r) = N{E^'^)e-'^^u\"^\^F, 



, , 1 i 

\m\ 1 , 1 + 2\m\Au 

II O / ' I I ' 

2 nj2||-E(2)| 



The potential being repulsive, there are no bound states in 3712 ■ 

The supersymmetry algebra now allows us to identify all the solutions of the supersymmetric 
spectral problem HipE = EipE from the eigenfunctions of Hq and H2 with non-zero eigenvalue. 
The key observation is that there are two kinds of non-zero (strictly positive energy) eigenfunc- 
tions ipE — Q'^i^E G Q^STi and ip^ = QipE £ QSTi because, if _E > 0: 

HQ^^pE = Q^H^E = EQ^E , HQijE = QH^e = EQ^e ■ 
The structure of the spectrum is as follows: 

• Ground states. 

There is a unique ground state -that belongs to STio and hence Bosonic- of zero energy 
4°^ = 0: ^S!ixi,X2) eKeiH. 

• There exist Q^-exact eigenstates of three types 

1. Q^-exact - henceforth, living in Q'^S'H - bound state eigen-spinors that belong to 

^2l)(^l'^2) = QHfi{Xi,X2) 
3 

2. Q^-exact - henceforth, living in Q^STi - scattering eigen-spinors that belong to STii: 

e['^ = E(°) > ^ , ^« ^ {x„x2) = gV2i) {xu X2) 

3. Q''^-exact - henceforth, living in Q^STi - scattering wave-functions that belong to 

sn2: 

i?(2) = E^^^ > ^ , ij% (xi, X2) = QH^^I {x^, X2) 

• There exist Q-exact eigenstates also of three types 

1. Q-exact bound states -henceforth belonging to QSTi- but living in ST-Lq: 

Ef = Ef^ , J > , i^fi{xi,X2) = Qi'^ (xi, X2) . 

2. Q-exact - henceforth, living in QSl-L - scattering wave- functions that belong to SI-Lq. 

i^(o) = ^W > ^ ^ ^S)(a;i,a;2) = Q^P^^,,,{x,, X2) 
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3. Q-exact - henceforth, hving in Q'^S'H - scattering eigen-spinors that belong to ST-ii. 



E^y = E^'^ > 



,(1) 



(2) 



h^ 



V^p(i)(a;i,X2) = QiIjU^^{xi,X2) 



Because the eigenfunctions form a total set in each sub-space we have the decomposition a la 
Hodge of the supersymmetric space of states: 




S'H = QS'H^Q^S'H^KeTH . 




Figure 1: The spectrum of the Kepler-Coulomb Hamiltonian (left panel). The spectrum of the 
supersymmetric Kepler-Coulomb Hamiltonian (right panel). 



3.6 Spin-statistics structure of the supersymmetric Kepler/Coulomb 
problem 

We have unveiled the spectrum of the planar supersymmetric Kepler/Coulomb problem solving 
the spectral problem of the two scalar Hamiltonians and using the supersymmetry algebra to 
obtain the eigenfunctions of the matrix operator. It is convenient, however, to look at the 
system as a whole, i.e., search directly for the spectrum of the 4 x 4-matrix supersymmetric 
Hamiltonian operator. 

With this goal in mind we define the "spin" operator 






h 



/ \ 
10 
0-100 

\ 0/ 



Clearly, [5*, A^] = 0, such that \S and N share eigenstates fulfilling a quantum mechanical 
spin-statistics theorem. The Bosonic eigenstates of N are zero spin eigenstates of \S, whereas 
the Fermionic eigenstates of A^ are one-half spin eigenstates of \S: 



-S 



( ^(0) \ 



V / 



-s 



( \ 






-s 



( 



\ 








\ 



h 

2 



/ 



V 








\ 
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Note also that neither the orbital angular momentum L nor the spin angular momentum 
S commute with H. The "total" angular momentum is the quantum invariant associated with 
simultaneous rotations of the Bosonic Xk and Fermionic il)k coordinateqj : 

J = L + S = -ih { xi- ^2^ j - ih ^^2 - ^l^ij • 

Therefore, besides the fact that [J,N] = 0, one can use [J,g] = [J,X] = to show that: 

[J,Q] = [J,Q^ = [J,H] = , 
where X is defined as in O): 

X = [g,g^]=I,-2N + 2g^g . 

We have two Clifford supersymmetric operators commuting with each other: H and J. The 
supersymmetric system as a whole is integrable. Now, the challenge is to find more Clifford 
differential operators commuting with the Hamiltonian. In [42j the authors found the super- 
symmetric version of the Runge-Lenz vector operator -henceforth, the supersymmetric KLPW 
vector operator-: 

One could guess the step from L to J and the need for the factor X is also no surprise given 
its role in the supersymmetric Hamiltonian H. A long computation ensures that the two 
components of this 4 x 4-matrix vector differential operator will indeed commute with the 
Hamiltonian and with the Fermi number operator: 

m, Q] = [Wk, Q^ = [Wk, H] = [Wk, N] = . 

Some work is also necessary to check that 

.. / . 2\ -^ 

[J, Wi] = ihW2 , [J, W2] = -ihWi , [Wi, W2] = -2ih [h~ — ]j 

Therefore, defining 

Ch = 



n? 



m-H) 



2 



the SO(3) Lie algebra is now closed -in the sub-space of states of energy in the range E G (0, ^^ 
by the Clifford operators Ci,C2,C^ = J 



[Ca, Cb] = -ihcabcCc , a, 6, c = 1, 2, 3 

and the Casimir operator is the 4 x 4-matrix differential operator: C^ = ^ 
the authors were able to find: 



cf + c| + a 



In 



'^^ 2^^^ h' [ (l-2iV)2 + 4CV ' '^t^ 2^^^ h^ [ i3-2Ny + AC\ 

such that Ej = ^ f 1 — .^.j|_^.^2 ) = Ej , the bound state energies paired through supersym- 
metry in STio and STii, reappear. 

^Recall that [ijXfc] = hskjXj and [S,ipk] = fefcjV'j, ^12 = —£21 = 1 , £11 = £22 = . 
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4 The planar quantum Euler/Coulomb problem and su- 
persymmetry 

Our second task is to build a supersymmetric quantum mechanical system inspired in the 
Euler/Coulomb problem: a massive/charged particle that moves on a plane under the influence 
of two flxed Newtonian/Coulombian centers, see next Figure. 




Figure 2: Location of the two centers and distances to the particle from the centers. 



4.1 The quantum Euler/Coulomb Hamiltonian 

The quantum Euler Hamiltonian is: 

/ - __^ A-— - — 
2171 ri r2 

Here, {xf = ±d, xf = 0) are the locations of the centers in the Xi-axis, and ai = a > a2 = 5a 
are the center strengths. With no loss of generality, we assume 6 G (0, 1] affording hetero-nuclear 
one-electron diatomic molecular ions. Thus, the strengths depend on the atomic numbers Zi 
and Z2 of the atoms and the 6 parameter is a positive rational number less than or equal to 
one: 



a = e Zi 



6 



Zi 



<i, , Zi,Z2eN* , z,>Z2 



Finally, the distances of the particle to the two centers are: ri = \/{xi — dy + x^ , r2 



Jun . 



Unlike in the Kepler/ Coulomb problem there is a parameter with dimensions of length in 
the system: the distance between the centers d. This allows us to use non dimensional spatial 
coordinates: 



Xi — )■ dxi , X2 — )■ dx2 , ri — )■ dri = dJ{xi — 1)^ + x^ , r2 — )■ dr2 = dy {xi + 1)^ + x^ . 

Note that there is also a fundamental action ymda built from the parameters of the system 
that provides a non dimensional Planck constant: h = -7^=. Assembling all this together, the 

V'mda 



linear momentum and Hamiltonian operators go to pi — ?■ 
new non dimensional operators are: 



-pi and Ji — > ^/i, where the 



Pi = -in-— 

OXi 



f n? ( 92 92 

/i = -— I 7^ + 



_1_ ^ 
2 \dx\ ' dx\) r\ r2 
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In this problem there is a non-obvious symmetry operator I2 — ?■ mdal2 where the non 
dimensional operator reads [16]: 



2 V / ri r2 

Just as the Runge-Lenz vector is quadratic in the momenta but unlike in the Kepler/Coulomb 
problem there are no more invariants in the Euler system which, accordingly, is only integrable. 
Explicitly, 

'Xi - IJTTT + 3^2 TT^ - 2X1X27^ 7; 3^17^ X2T. H h 



dxl dx\ 9X19X2 dXi (9X2/ ^1 T2 

and a little algebra shows that: [Ji, 12] = /1/2 — I2I1 = . 

4.2 Separability of the Schrodinger equation in elliptic coordinates 

Because of the quadratic in the momenta symmetry operator, we expect that the Schrodinger 
equation will be separable in some coordinate system on the plane. To skip the singularities in 
the centers one can cover the plane by two open charts: the first chart is the open set M^ x N 
in ]R^/{(— 1,0)} = M"*" X S^ where N is the open North hemisphere in §^. The second chart is 
the open set M"*" x § in M^/{(1,0)} = M"*" x S-*^ where § is the open South hemisphere. Both 
charts must be glued at the abscissa axis X2 = 0. Totally adapted to this topological situation 
are the elliptic coordinates; the half-sum and half-difference of the distances to the centers of 
the particle: 

u = -{ri + r2)^{l,+oo) , ■(; = -(r2-ri) G (-1,1) 

that parametrize a two dimensional infinite strip ¥?. The Cartesian coordinates are obtained 
through the change 



Xi 



uv G ( — 00, +00) , X2 = ±a/(m^ — 1)(1 — f^) G ( — 00, +00) 



which is a two-to-one map -one per chart- from M^ to ¥? except at the X2-axis, which is one- 
to-one mapped at the boundary of E^: d¥? = {{u = l,v), {u,v = ±1)}. 
The Euler Hamiltonian in elliptic coordinates is of the separable form 

h = -1. (Hu + Hi 



and the symmetry operator also separates 

i 1 



y2 _ y2 



u^ - l)H, - (1 - v'')H^ 



The ansatz iJje{u,v) = riE{u)S,E{v) converts the spectral problem 

iiipE{u, v) = E^e{u, v) (6) 

into separable: 

-¥{u^ - 1)^(m) - n\^{u) - [2(1 + S)u + 2u'E] r]Eiu) = 1^^) (7) 

-¥{1 - v')^{v) + ^'^^(^) + [-2(1 - S)v + '^v'E] Uv) = -lUv) • (8) 
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The Schodinger PDE equation (|6| becomes the two coupled ODE's ((7|-([8]) where the separation 
constant / is the eigenvalue of the symmetry operator / = — 2/i — 2/2. 

We could try to solve ((7])-([8]) directly but we still perform the following change of variables: 

1 1 r ^T 

X = -arccoshw G [0, 00) , y = -arccosf G [0, —J . (9) 

Equation ([T]) becomes the Razavy equation (10) [H], and ([s]) becomes the Razavy trigonometric 



(11) or Whittaker-Hill equation 

+ {C cosh2x - MY r]Eix) = X^six) (10) 

+ {/3cos2y-Nf^E{y) = f^Uv) ■ (H) 



(PrjEJx) , ,^ „„„,, ^^, ,^^2 

dx'^ 



dy"^ 



The parameters in the Razavy equations (10) and (11) are defined in terms of the energy and 



the eigenvalue of the symmetry operator in the form: 



2^/-2E 


h 

2^/-2E 



We stress the following subtle point: the Razavy equations are defined for fixed C, M, and A 
or for fixed /3, N , and /x. We obtain, however, Razavy equations for parameters determined 
from E and /. Therefore, we address an infinite number of entangled Razavy and Razavy 
trigonometric equations. 

In Reference [38] it is shown that the Razavy and Whittaker-Hill equations for M and 
N positive integers are quasi-exactly solvable (QES) systems and all the finite solutions - 
polynomials times fast decreasing exponentials- are found by algebraic means. Our strategy 
will be to use this information in the search for the bound state eigenvalues and eigenfunctions 
of the Euler/Pauli Hamiltonian. Our results have been partially published in [52] . Thus, we 
shall not repeat the analysis here. Instead, we shall develop the program in the supersymmetric 
version of the Euler problem. 

4.3 The supersymmetric modified Euler/Coulomb Hamiltonian 

In [52] we gave arguments for selecting the following superpotential 

W{xi,X2) = -^{ri + 5r2) (12) 

in order to develop a supersymmetric quantum mechanical system from two fixed centers con- 
taining a mild deformation of the Euler/Coulomb system in the A^ = sector. In fact, from 
the superpotential partial derivatives 



dxi 2h \ Ti r2 7 ' 0x2 2h \ri r2 

d'^W _ 1 f/1 {xi-lf\^^ ,^ /[ ^ (xi^p2 

,3 



dxl 2h 1 Vn Arl J'^ ^^^ [r2 Ar? H^ + '^^ 



^ = -^{u-||)^'-^) + U~l^'^'+^) 



2, 

2t 



aw _ 1 / X2(X1-1 X2(X1 + 1) \ _ d^W 

9xi9a;2 2h \ Arf 4rf J dx2dx\ 
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we obtain first the supercharges. In turn, the scalar Hamiltonians, 



H. 



Ho 



--A 
2 



2 



W , 2 





4 Y 


1 5 


Vr2 ri 


4 Y 


1 5 

+ — + — 



and the matrix Hamihonian: 



Hi 



Ho + H 



-(xi+l)^-xl 



3:^2(3:1-1) 



+ 6 



3:2(3:1+1) 



3:2(^1-1) 



^5— +S^-zs 



3:2(3:1+1) 

r: 



m 



H2] + 



(3:i-l)'-3:| _^ ^(xi+l)2-x2 



'2 



are derived. Now, the rationale for the choice of the superpotential (12) is clearer: at the limit 



where the two centers are superposed, d = and ri = r2, the superpotential, the supercharges 
and the superHamiltonian become those of the supersymmetric Kepler/Coulomb problem (with 
non-dimensional strength 1 + 5 instead of 1). 

In this one-parametric deformation of the Kepler problem a very subtle conundrum arises. In 
the Kepler/ Coulomb case, the non-supersymmetric K and the A^ = scalar Hq Hamiltonians 
only differ in the shift by a constant necessary to push the ground state energy to zero as 
required by supersymmetry. In the Euler/Coulomb case, /i and ^0 differ in a non-constant 
potential such that in the Kepler limit becomes the constant shiftQ 



Vs{xi,X2) 



~1^ 



ri , r2 



4t/2 



l + 6' + 6[- + -- 

^2 ri rir2 



lim Vs{xi,X2) 

d-s>0 






The role of this potential is to shift the negative bound state energies in a harmless way: Ho, like 
Ji, still admits a symmetry operator that is quadratic in the momenta and the supersymmetric 
spectral problem in SHq is still separable in elliptic coordinates !! In other words, we choose the 
superpotential in such a way that the separability in elliptic coordinates of Hq is preserved even 
though we must add a "classical" piece -important when h tends to - to the Euler/Coulomb 
non-supersymmetric Hamiltonian. 

In fact, Hq and the symmetry operator I2 in elliptic coordinates are still of the form 



Ho 



1 



u^ 



^1°' + ^1°) 



m 



u^ 



M 



m^:' 



^)^r 



where Hu and Hi, are now 



B 



u 



5! 
2 

2 



du^ du 

/I _ 'i\_ 

dv"^ dv 



+ j-,{l + 5)\u- 



{l + 6)u 



+ j-,{l-5ni-v')-{l-5)v 



(13) 
(14) 



The separation ansatz ^^ (M,t>) = ^e vA^e \") plugged into the supersymmetric spectral 
problem in the A^ = Bosonic subspace 



.,(0) 



Hoi^';\u,v) = Ei;'^\u,v) 



(0), 



^We temporarily come back to dimensional coordinates in order to see the limit. 



(2) 



Obviously, the same situation happens in the N = 2 sector. H2 and /j are given in the same way in terms 
of Hu and Hy that differ from Hu and Hij in the sign of the last terms. 
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reduces the PDE Schrodinger equation to the system of separated ODE's 

2 ^ / n I X\2 









d°^(^^) 



Irifiu) (15) 
-/d°^(^) (16) 



where the separation constant I is the eigenvalue of the symmetry operator / 



2Hq - 21. 



(0) 



4.4 Bound states from entangled Razavy and Whittaker-Hill equa- 
tions 



The change of coordinates ^ transforms ( 15 ) and ( 16 ) respectively in the Razavy and Whittaker- 
Hill (three-term Hill) equations 



72^(0) 



"^ !f i""^ + (C cosh 2a: - Mfr^f{x) = \r^f{x) 

^^ + (/3cos2y-iV)^4^y) = ^^if{y) 



where the parameters are now: 



(17) 
(18) 



^4vl(^+^)^-^^ 



^-hV 1(^-^)^-2^ 



M- _ 2(1 + 5f 4 4(1 + ^)^ 

^^2 2{l- 5f ,,2 ^^r 4(1-5)2, 



4.5 Eigenfunctions from the Razavy equations 



In [38] it is shown that the Razavy equation (17) is a quasi-exactly solvable algebraic equation 



that admits n + 1 finite (polynomial times decaying exponential) solutions if M = n + 1 and 
n e M is a natural number. Moreover, there are ra + 1 solutions for ra + 1 different values of 
A characterized by an integer m = l,2,3,---,n + l. All the eigenvalues between and -^^ — — 



n? 



of the supersymmetric modified Euler/Coulomb spectral problem in SI-Lq are obtained in this 
way: 



--'- - -'^^'-Vi^-jM 



Concerning the eigenfunctions, we search for solutions of the M = n + 1 Razavy equations by 
means of the series expansion Q 



/ N " 'J"l'r4-l'l V-^ ( — l)'^-Pfc(A) 



,2x 



k=0 



{2C)''k\ 



Cn 



4(1 + 5) 



where Pfe(A) are polynomials of order /c in A to be fixed. The ODE Razavy equation is then 
solved if the following three-term recurrence relations among the polynomials hold: 

Pfc+i(A) = (A - {Ak{n -k) + 2n + l + C)) Pk{X) - (4A;(n + 1 - k)C) Pk-i{X) , k>0 . 



^This ansatz, and others related to this, aUows to represent the Razavy Hamiltonians in terms of differential 
operators that belong to the enveloping algebra of §L(2,M), see [38] and [37] , 
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In particular, if A„m, m = 1, 2, ■ ■ ■ , n + 1, is one of the n + 1 roots of Pn+i, Pn+i{^nm) = 0, then 

and the series truncates. The degeneracy in the energy is broken by the eigenvalues / of the 
symmetry operator I2 provided by the roots \nm in the form 



h^ 



(Anm - {n + if) - 



4(1 + ^)- 

^2 



which distinguishes between the different polynomials Pmi^nm)- 
We solve the finite-step recurrences for the lower-energy cases 



n = 0: 



Co = ^^^, Po(A) = l, Pi(A) = A-(l + 



A, 



01 



1 + Co^ 



n = 1: 



Ci 



2(1 + ^) 



Po(A) = l , Pi(A) = A-(3 + C?; 



P,(\) = A2_2(3 + Ci')A + 9 + 2Ci' + Cf 



An = 3-2C1 + C1 



A12 = 3 + 2Ci + Ci 



n = 2 



C2 

^2(A) 



A 
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4(1 + ^) 



^o(A) 



Pi(A) = A-(5 + C2') 



(A-(5 + C2'))(A-(9 + C2^))-8C2^ 

(A - (5 + CI)) [(A - (5 + C|))(A - (9 + CI)) - 16CI] 

A22 



Cl + 5 



CI 



7-2./1 + 4C 



A 
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CI 



7 + 2./1 + 4C 



and show the results for the lower eigenvalues and eigenfunctions in the next Table: 



Energy 



Eigen-function 



Separation constant 



^0 



n'i'^\u)-. 



^(0) ^ 3(1 + 5)^ 
^1 2S2 



(1 + .^)" 



3(1 + ^)" 



-(l + «) 



/^ ' — — il— 



^^<a#i^ + (1 + ^) 



j(0) ^ 16(1 + 5)^ 
^2 9R^ 



''22 '."•' — 4(l + a) 



, 2(i+a) 



'723 (") ^ 



V^u^ - 1 
3R^ V OR* 



3fi'^ ~ 3k2 8(1 + ,5) / 64(l + a)2 



r(0) 



= 2 32(1 + 5)^ 



-4^-51^J^- 1^/64(1 + 5)^ +9M 
-¥-^^^^ + ^v'64(l + .)2+9ft^ 
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4.5.1 Contribution from the Whittaker-Hill equations 

For these values of the eigenvalues En and Inm two of the parameters of the Whittaker-Hill 



equations (18) become: 






but the other one Nn = , "^ = is not a positive integer if n > 1 and 5 = ^. The 

existence of solutions of the equation 



Zi 



{k + 1) v/n(n + 2) - (n + l)^/k{k + 2)1 = Z2 [(A; + 1)^^2(72 + 2) + (n + 1) v/A;(A; + 2) 



over positive integers n and k would be necessary to simultaneously find M = n + 1 and 
N = k + 1. It is clear that this is not the case and there are no finite solutions other than the 
ground state in the supersymmetric two-center problem. The analogous equation in the non 

supersymmetric case is: 

_ n+k+2 

n — k 
which admits solutions for positive integers n and k found by Demkov and his colleagues over 
forty years ago, see [51] and, e.g., [55] . 

li n = 0, Eq = 0, /q^ = 0, however, the Whittaker-Hill equation is also QES, see [39], with 
a unique finite wave function: 

The analytic wave function of the ground state in the SHo sector is: 

£f=0 , .<)(„. ,,)=,S'MC(")=exp[-?<ii^l.exp[Hii^l ; (19) 

thus, there is a normalizable Bosonic zero energy ground state, a zero mode ipQ^ {u, v), in the 
supersymmetric two-center system. Supersymmetry is not spontaneously broken. 

In the WH equation for the above parameters we could try a solution of the form [39j : 






w'' , w = e^'y 



fc=0 



which solves (18) if the "recurrence" relations between the Q-polynomials hold: 



Qfc+i(/i„) = (^„ - {4k{Nn - 1 - fc) + 2iV„ - 1 + (3l)) Qk{iin)-{4k{Nn-k)(3l) Qt-iM ,k>0 

This strategy, however, is not useful in this situation because the WH equations are not QES 
if n> 1 (A^^ A; + l). 

Instead, we consider the WH equations in their algebraic form: 

dv"" ~~ l-v^ dv l-v^ ^""^^^ ~ 

In^ A{l-5f _ 2(l-(5) __l/p 2(1-5)^ 
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Now, following the standard theory of Hill equations, see [57j, we search for power series solu- 
tions: if A; G N is a natural number. 



eiiv) = E 






.(2) 



A. 



nm JO) 



.(3) 



1 



nm ? ^nm 



\\^nm ^J'^nm ^ '^nm'^nmj 



fc=0 



.(fc) 



k{k 



^{(A_-(fc-2)2)cit'^ + 5_cit')+C_cit'^} , fc>4 (2Q) 



We encounter fourth-term recurrence relations, a difficult situation to deal with, although the 
two basic solutions are easily characterized: 






.(0) ^ 2 r^^^ = 0- f"-"' iv\ = 1 + V°° r^ 

■'nm+ ^1 ^nm+ "• Snm+V^/ -^ ' /—dk=2 ^nm+ 



2. a 



(0) 



c^'^ 



,{0) 



oo (fc) 



1: ^nm-(^) =^ + Er=2 C. 



V . 



These series converge in the open (—1, 1) interval and are extended to cover the singularities 
V = ±1 setting the values: 



d°^+(i) = i + E^^ 



(fc) 

nm+ 



kjk) 
^"nm-\ 



k=2 

oo 



d°l(i) = i + E^^ 



(fc) 

nm- 



k=2 

oo 

Snm— \ ^) -*- ~r / j \ ^) ^nm 



(21) 



k=2 



k=2 



The series for the ground state, for instance, are easy to find: 



.(0) 



^01 (^) = exp 



2(1-5), 



dv 



2(1-5) 



^2 



-exp 



fc=0 

2(1-5) , 

^2 






-01 



1 /^2(l-5) 



A;! 



^2 



5 "-^01 



ci^^=eff(0) = l 



fc=l 



(fc)^ fc-1 (1) _ '^6 



(0) 
01 



01 



-01 



dv 



(0) 



2(1-5) 

^2 



eJ?(±i) = E 



(±1)^ /2(l-(5) 



A;=0 



k\ 



n? 



exp 



2(1-5) 



Any other solution of the WH equations is obtained by specific linear combinations of the two 
basic solutions. We will choose linear combinations of the general form 

ej.(^) = c^i^lM + C-d°L(^) (22) 



In fact, any choice of {c^,cJ) G £? in (22) fixes the extension to the boundary of the elliptic 
strip E = [1, +oo) X [—1, 1] of the wave functions il)nm{u^v) = ^nm{u)^nm{v) . Id est, extensions 



(0) 



',(0) 



of Hq in L (E) are determined from the values of ijjnmiu^v) at the boundary: ^nm(l,f) and 
'4'nm{u, ±1). We remark that these extensions are not essentially self- adjoint in general; different 
eigenfunctions have a small overlap for generic values of h. 



4.6 Two centers of the same strength 

In the case 6 = 1 when the strength of the centers is identical the supersymmetric spectral 
problem in SHo simplifies remarkably: 



-¥{u' - 1):^ - h'u4- + f^(^' - 1) - 4« - 2Eu' 
dw^ du V II 



.(0) 



% H = IVe (u 



(0), 



-h\l - v^)$^ + Wv^ + 2Ev^ 
dv^ dv 



,(0) 



(0), 



eE\v) = -leE'iv 



(23) 
(24) 
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Again, (23) is equivalent to a Razavy equation 

(fr]f{x 



2 JO)/ 



dx'^ 



+ (Ccosh2x - MYt]'-^\x) = X r]'-^\x 



S^)i 



(25) 



with parameters 



c 



16 



-2E 



M^ 



after the change of coordinates: u = cosh2a;. 



X = M^ + ^{I + 



h? 



16, 



The v-equations (24), however, become the Mathieu equations 



(0), 



dy' 



+ (acos% + a)eF(y) = 



(26) 



(j= TTj(/ + ^) 



with parameters 

_ AE 

under the change: v = cos2y. The strategy to solve these two entangled equations is the same 
as in the case of two centers of different strength. 

First, we search for finite solutions of the Razavy equation. The procedure is identical and 
we only need to replace (5 by 1 in the formulae of the previous Section §. 4.5. We now have 



E(°) = ^ I 1 



1 



in + r 



Cn 



h' 



h^n + l) 



{Xnm-{n + lf) 
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where A^ 



m = l,2. 



,n + l are the roots of the polynomials P„+i(A) that cut the series. We 



thus show the results for the lower eigenvalues and eigenfunctions in the next Table: 



Energy 


Eigen-function 


Separation constant 


Ei°) = 


^S'M = e-S 


7-{o) _ n 
-'oi — u 


p(0) _ 6 

^1 - w 






( r(0) h^ 12 r, 
^11 --X" ft2 -^ 

< 

7-(0) ri^ 12 I 9 

I -'12 - ~x F" + ^ 


< 


[ rjt,'{u) = -e-Vi^2{u-l) 


p(0) 64 
^2 9fi^ 


. 


f 'n^ ■i" 


< 


f 7-(0) _ fc2 128 
-'21 " 9^2 


ri^^l>{u) - -2e 3h^ y/u^ - 1 


^ 4-11. 


^^2^ = -f-il-iV256 + 9fi4 


[ n\l>(u)=^-^e SS^^_1^„_1_^1+|||J 


i /^3^=-f-# + iV256 + 9;.4 



Nothing new with respect to the non equal centers case. 

4.6.1 Contribution from the Mathieu equations 

The novelty comes from the Mathieu equations: unlike the Whittaker-Hill equations these 
equations are never quasi-exactly solvable -there are no finite solutions whatsoever- but, instead 
there are lots of solutions that can be described analytically in terms of the Mathieu sine and 
cosine special functions. 
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As in the non-equal centers system, the ground state is an exception. It is not ruled by any 

p-E'o = 0, o"oi = p( 



Mathieu equation, oq = p-^o = 0, o"oi = A(/oi + Eq) = implies that: 






dy' 

If the two centers have the same strength an important discrete symmetry arises: the ri -H- r2 
exchange is not detectable, v O —v, or, y f-)- y + | is a symmetry of the system and we remain 

with the only invariant solutions Qi (y) = B under this transformation. There are only two 
independent choices: B = 1 and 5 = 0. The first choice is "even" in f , the second choice 
is "odd" in v but negligible. The zero energy ground state is therefore built from the even 
^-independent wave function: t/'q^ (M,t>) = ^i ('w)^oi ('^) = ^~^ ■ 

The parameters of the Mathieu equations determined by the spectral problem for positive 
energy {n > 1) are: 



f9 ' '"" h? 



The V -H- —V symmetry of equation (24) is translated into the y ■^ y + |n, n E "L, infinite 
discrete symmetry of the Mathieu equation. The Mathieu cosine and Mathieu sine special 
functions are obtained from the Bloch-type solutions of the Mathieu equation ruled by the 
discrete translational symmetry with Floquet indices determined from the parameters a and a. 
Because the Mathieu equation is blind to the v ^ —v exchange, if ^nmiv) is a solution ^nrn{—v) 



also solves (24). In [52] we chose even and odd in v combinations of the Mathieu functions, a 
situation closely related to the hidden quantum supersymmetry in Bosonic systems unveiled in 
j51j . Instead, here we choose the combination 

^im(^) = C [anm, Qn, arccos(t;)] + iS [anm, qn, arccos(t;)] , (27) 

to build positive energy eigenf unctions, a choice designed to go to the Kepler/ Coulomb system 
at the d = limit. Here, C[a,q,z] and S[a,q,z] are respectively the cosine and sine Mathieu 
functions and (?„ = ^, anm = — ^^. 



We now list some of the lower values of the parameters 

A + i- + 1 
a -'^^1 

«12 - 14 - I2 + 



^i = F ' ^!":f!ft! 



rj4 rj2 I 4 



12 - 9f,4 , S "21 - J- -^ 9^4 , 



"23 = ^ + ^-^^256 + 9^4 

and offer a Table showing the probability densities of some eigenfunctions for several values 
of h. For instance, h = 0.7 is the value of this non-dimensional parameter for the hydrogen 
molecule ion. h = 114,7 corresponds to the ionized helium atom; note that the radius of the 
nuclei is of the order of 10~^^cm, etcetera. One notices that the smaller h is, the more classical 
is the system, the wave functions being more concentrated around the centers. 

4.7 Comparison between the J\f = 2 super symmetric and the A^ = 
spectra 

The energy spectra of the Euler Hamiltonian - given by the solution of the system of equations 
(71), (tel)- and the Hq Euler supersymmetric Hamiltonian reduced to the SHq sector of the 



supersymmetric Hilbert space -characterized by the system (15), (16)- are listed below 



X = , A/' = ; 

n^A-^<1 F - 2(1+^)^ p(0) _ 2(l+,5)2 , . 



Zi — -^ ' -^" h^{n+iy ' -^^ n? y- (n+l)2 
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Table 3: Af = 2 supersymmetric wave functions in SHq of two equal centers: S = 1 



\i^^l{xi,X2)\^ h = 0.7 



n=0, m=l 



n=l, m=l 



h= 1 



h = 2 



h = A 



't--: 








n=l, m=2 




n=2, m=l 




I'.-n 





7 



n=2, m=2 




.-^v-^^,. 



^:;? 





n=2, in=3 
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The degeneracy of the nth level is n + 1 in both cases and it is split by the eigenvalues of the 
symmetry operator, respectively / and I^^\ 

As in the supersymmetric Kepler/ Coulomb problem there are bound states in SHi of the 
form: 

^i'lin,r2) = Q^tP^iin,r2) , n>l, m = 1,2, ■ ■ ■ ,n + 1 . 

Unlike in the supersymmetric Kepler problem, there is one zero energy wave function in STii. 
Coming back to dimension-full coordinates and parameters, it is: 






gVo'^(xi,X2)=0 

/ 



E, 



(1) 







V 







£2 
4c(ri 



e A^ 



(28) 



/ 



here, the " — " sign occurs for 0:2 > whereas the " + " sign arises when X2 < 0. In the limit 
where the two centers coincide this wave function becomes: 



/ 



lim'?/'A^^(a;i,X2) 



\ 

X2 



XI 







e-5=^f±^v/^ 



and we observe that its norm diverges: 



lim 

d-s>0 



dxi dx2 



—00 J —00 



,(1) 



^0 (2;i>a;2) 



lim27riro ( ^^ (l + i^) 



d-s>0 



n? 



ij 



Amad 



V h^ 



1-6) 



+ CXD 



where Kq{z) and Io{z) are respectively the first-order modified Bessel functions, see |J^. Thus, 
we confirm that there is no fermionic zero mode in the SUSY Kepler problem. Note that ground 
states in that problem must live in the scalar representation of SO(3). 

The comparison between the Bosonic and Fermionic zero modes can be seen in Table 4, 
where the center on the right is twice as strong as the center on the left and 5 = |. In the 
Bosonic ground state the electron is concentrated around the stronger center for small h but 
becomes spread over the two centers when h increases. In the Fermionic ground state the 
superparticle behaves in the opposite way!: for small h it is concentrated in the weaker center 
and two peaks on the two centers arise for larger h. Moreover, for any value of h the probability 
density of the ground state in STii is always peaked at the centers. In any other respect, i.e., 
concerning the scattering solutions with energy greater than p(l + 5)^, the structure of the 
spectrum is qualitatively identical to the spectrum of the supersymmetric Kepler problem. 
There are scattering states in STio paired via the Q^ supercharge to scattering states in STii 
and scattering states in 57^2 paired via the Q supercharge to scattering states in STii. AH this 
is depicted schematically in Figure 3: 



5 Two centers collapse in one center 

In this Section we shall analyze how the two main spectral problems described respectively 
in sections §.3 and §.4 are connected. The link appears at the d = limit of the two-center 
problem. In order to go to this singular limit of the modified Euler/Coulomb problem it is 
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Table 4: A/" = 2 supersymmetric zero modes. 



Probability density /i = 0.7 



\i>'£'{xi,X2)\^ 




h=i 



h = 2 





h = A 





5=1 




E 


H 





|| 












2(1+5/ 









«o ^ 


. Hi , 


, Hn 








E 

2(1+5)'' 


i 




11 


■r = 








h 
































n 









Figure 3: The spectrum of the Euler Hamiltonian (left panel). The spectrum of the supersym- 
metric Euler Hamiltonian (right panel). 



29 



necessary to restore full dimensional variables. Thus, with no change of notation the parameters 
and physical variables to be dealt with in this Section have the proper dimensions. We shall 
perform the limit in the case 6 = 1 because we have full analytical information for two equal 
centers. 

• We start from the ground state characterized by: n = and m = 1. The eigenvalues of 
the Hamiltonian and the symmetry operator, the Mathieu parameters, and the ground 
state wave function, are respectively: 

At the (i — )■ limit we have: 



r = lim-u = \/a;? + Xn , cp = limarccos(-) = arctanf — ' 
Therefore 



-') 



',(0)/^ ,„X _!• „,,(0). ^ ima\ ^2rna. 



and the ground state of the one-center problem with twice the electric charge appears in 
the (i = limit. 

We next consider the doublet labeled by n = 1, m = 1 and m = 2. The eigenvalues of 
the Hamiltonian and symmetry operators are: 



^r 






^' ' Iu = -'i + 2{mda)-'-^^ 



The corresponding Mathieu parameters and eigenfunctions read: 



Qi 



3{mdaf an - 4 + ^^^ + -^^ 



^4 ' ^ 

4 '^~fp ^ ¥ 



1 D (mda) I 6{mda) 

^12 — 7 — ^— E2 1 ki 



"max 3/2 _2mc«. 



tl^f^iu.v) = i — \ e h2 « ^2{u + (i) |C[aii, gi, arccos(-) + iS[aii, gi, arccos(-)] 



^12 (W)^) = - (-p-j e "^?^" a/2(m - d) |C[ai2, gi, arccos(-) + iS[ai2, qi, arccos(-)] i 

Because 

1 ff 

hmgi=0 , aifc = lim aifc = - , A; = 1,2 , lim/n = lim/12 = --- 

d-fO d-fO 4 d-5>0 d->0 4 

V V 

cosy/aik'^ = limC[aifc, gi, arccos(-)] , siuA/aifc ip = hm S[aik, gi, arccos(-)] , 

d-^O d d^O d 



we find 



/ ~ \ 3/2 
,(0)/ N ,. ,(0)/ X (ma\ ^2 .m&r ilu, 

2 2 d^o \ n^ I 






,(0) / N ,. ,{0)/ N ma\ W2 ^m^r -i\u, 

ip\\[r,ipj = hm.ip{2 [u,vj = — I -—- r ' e ^^ e ^^ 



Again we find that the eigenvalue of the Hamiltonian operator remains the same at the 
(i = limit, the eigenvalues of the symmetry operator go to the eigenvalues of the square 
of the angular momentum, and the wave functions become the eigen-functions of the 

Kepler/ Coulomb problem with spin of one-half. Note that we have used: \ \ = ±7^ 



4 2' 
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Finally, we consider the triplet n = 2,'m = l,'m = 2 and m = 3. The eigenvalues of the 
Hamiltonian and the symmetry operator are: 



t-2 _ I28(mda)^ 



E. 



(0) 



9h^ 



'22 



'23 



hl-- 

h2 128(m.da)^ ^2 



" 2 



9fi2 
128{mdaf 



256(mdQ)2 



1 + 



9^2 



I S2 /-I I 256{mda)2 

+ TV^ + ^ 



whereas the Mathieu parameters read: 



Q2 



32{mdaf 
9hF~ 



^22 



a2i = H 

1 1^ 64(mda)2 



64(mdQ)^ 
9F 



+ 



9h'^ 



256(mda)2 
9ft4 



'^23 = 

The eigenfunctions are more complicated 



1 _|_ 64(mda)^ _ i / -, 

2 "•" 9/1* 2 "^ ' 



256(mda)2 
9F 



(0), , /mQ;\2 _4ma„ ^^r -r f f t) 

"^21 (""? "^J = I "^ j e ^h^ VM^ - a^ <^C[a2i,5'2,arccos(-) +26[a2i,g2,arccos(-)J 



^J?(M,t;) 



- e 'm^ 



/maV- 

y'W) 8 



mo; 



IQma 



-M- 1 + 



256(m(ia)" 



',(0) 



^23(m,w) 



3h^ V 9^2 

X < C[a22, g2, arccos(-) + iS'[a22, ^2, arccos(-)] 
ma\2 3h^ 



/ma\ 

\ h? J 8ma 



4,ma 

e ^ft^ 



IQma 



-u-1- 



25Q{mda)'^ 



X < C[a23, q2, arccos(-) + iS'[a23, q2, arccos(-)] > 



but the limits are similar: 
lim q2 = , 
lim /21 = lim /22 

d^O d^O 



0:21 = lim a2i = 1 = hm 022 = 022 

d-5>0 d-5>0 



0:23 = lim 023 = 

d-5>0 



-n^ 



lim I23 



, A; = 1,2,3 



cosa/q^V? = limC[a2fe,g2,arccos(-)] 



sm^/a^(p = limS'[a2A:,g2,arccos(-)] . 
d^o a 



Thus, we find 



,(0)^ 



',(0)^ 



^11 (^,V^) = limV^2i (w,^) 

a— >0 



~N 2 

ma \ _ 2mc 

_J re ^ 



3*'/' 



,(0) 



■,(0)^ 






ma 



r e~l^'' e-'^ 



•SO), 



'M, 



Wwir,(p) = lini^y(M,w) 

a— >3 



m,a\ f 3h^ ^ \ _2ma. 
— - — - - 4r 1 e ^^ 



again falling in wave functions of one doubly charged center, in this case with angular 
momentum 1. 

3D graphics of this analysis are shown in Table 5. We remark that the graphics are drawn in 
non-dimensional variables and some of them cannot skip the 1 in the d = limit displaying 
non-smooth tendency to the Kepler wave functions; namely the n = 2, m = 1 and n = 2, m = 2 
wave functions. 
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Table 5: 3D Plots of the probability densities: Kepler versus Euler. 



Probability 
density 



h = i 



j = 0, m = j = ^, m = ±^ j = 1, m = ±1 j = 1, m = 



Kepler 





^=*^ 




h» 1 



n = 0, m = 1 n = 1, m = 1 n = 2, m = 1 n = 2, m = 3 



Two Centers 






h» 1 



n = 1, m = 2 n = 2, m = 2 



Two Centers 
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5.1 Collapse of two centers of different strength 

When the two centers have different strengths, 5 7^ 1, the one center collapse happens exactly in 
the same way as compared to the equal two-center collapse because lim^^^o Bnm = lini5=i Bnm = 
0. Id est, when the charges are superposed the only things that matter is the total charge. We 
now offer this analysis for completeness which will shed light on the physical nature of the 



recurrence relations (20) and the basic solutions (21) of the Whittaker-Hill equations. Because 
we prefer to deal with power series with purely numerical coefficients, we shall stick to non- 
dimensional f-variables and put all the physical dimensions back in the w-dependent part of 
the wave functions. We recall that: 

A --^(t I 4(1 -^)^ \ r. _ 2(1-^) __^/ 2(1-5) 



are the non-dimensional parameters of the algebraic WH equations. 

• Start with the ground state: n = 0, r/i = 1. The values Aqi = — f f^ ) , -B, 



01 

(fc) 



^^2 , and Cqi = ^ ^4 ' means that the recurrence relations (20) are solved by Cg^ — 

^ i p ' ] . Thus, ^Qi (f ) = ^^0 h. ( P I ^^ must be multiplied by the unique finite 
solution of the Eq = Jqi = Razavy equation: 

,(0)/ ^ f'ma{l+6)\ 2ma{l + S) ^ 2niad{l-S) 



%i'iu,v) 



h^ 



It is clear, like above, that now 



(0)/ ^ _ I f^OL \ _2ma ^ _ (0), 



and we find the ground state of only one-center with zero energy and momentum and 
electric charge {1 + 5)a. 

[ 772 ^ 1 

Next, we consider the doublet: ra = 1 , < „ . The eigenfunctions are of the form 

y m = 2 



h^ 



3/2 



where the ^ji (f ) and Ci'ii'^) ^"^^ linear combinations of the basic solutions of the WH 
equations to be identified in such a way that the Kepler/Coulomb wave functions of 
energy Ei = ^^° and angular momentum | and — | arise at c? = 0. 

The pass to the limit is identical to the collapse of two equal centers in the w-part of 
the wave function. Therefore, we will describe in detail only the collapse of the v-part. 
Because the d = limits of the WH parameters are 



1 



lim All = lim A12 = - , lim Bn = lim B12 = = lim Cu = lim C12 



the recurrence relations become very easy: if we denote q = 1,2, 



1 (n\ rs^i 1/1 

"6 U 7^'" ' ^"^ k{k-l) V4 



(2) _ __ (0) (3) _ _-L / 1 _ -, \ (0) (fc) _ -^ / -^ ^,^ r,\2 \ {k-2) 

^19- 0% ' %- «U M ^19 ' %~ 1.(1. -^A A ^'^ ^> r^i9 
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The basic solutions, corresponding respectively to the choices ci^j^ = l,c^g_^ = and 



Cig_ — 0) Cig_ — 1, are: 



J2) _ _^ (4) 0_ (6) _ ^^ ___ (0) , . _ >r- (2k) 2k 

'^lq+ g ' '^lq+ ]^28 ' ''^ 1024' ' ^"^1+^ ' Z-^ ^^(1+ 

k=0 
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.(3) _ i „(5) _ ^ (7) _ J30_ (0) , , ^ -^ (2k 

-iq- g ? ^iq- ^2g ' i9- 1024' ' ^1?-^ '^ Z^ ^l'^^ 

fc=0 



We know that C[l/4, 0, arccosw] ±iS'[l/4, 0, arccosw] goes to exp[±i|] at c? = 0. Compar- 
ison with the power series expansion 

C[l/4, 0, arccosf ] + iS'[l/4, 0, arccosf ] = 



V2 V2 V2 V2 V2 



\256 256/ Vl024 ^ 1024/ ^ 



V2 V2 

and the analogous series with the relative minus sign, tells us that the right combinations 

;|] and lim^^od" (^) = exp[-?|] 



such that limrf^cCii ('^) = ^^Pi^f] ^^^ lii^d->o ■^12 ('^) = ^^pI^"^!] ^^^- 



?lg l^J - /^ ^lq+ ^ 2^2 ^^~ ' 

Henceforth, 

/ ~ \ 3/2 
lim<(«,t;) = limr/S?(«)ei?(t;) = ( ^ ) r^^ g-^^ e^l^ = ^^^(r,^) 



, A ^11 V"^' ^/ V'„'/li V^/Sii v^/ 1 »-2 / 



2 2 



3/2 



lim^i?K^) = lim^i?K^)d?M = -(^) r^/'e-l^^e-*i^ = ^(°\(r,v9) . 

d-5>0 d-5>0 \ h'' J 2~2 

m = 1 
Finally, we address the triplet state n = 1 , ■^ m = 2 . The c? = limits of the WH 

^ m = 3 
parameters are in this case 

lim A21 = hm A22 = 1 , hm A23 = , hm B2„ = = lim C2iq , g = 1, 2, 3 . 
The basic solutions for g = 1,2 are: 



16 ^'^^ 128 



c£-=0, 4;L = 1, 4^=0 ,VA;>2,--- 



oo 



fc=0 fc=0 

Comparison with the power series expansion 



if^ iv'^ iv^ 5iv^ 7w^^ , iix 

C[l, 0, arccost;] + iS\l, 0, arccoswj =^ + 1;-^ _____ __ _ ___ j^ Q [y ) 

show us that the combinations 

d? (^) = ^diV (^) + d?- (^) , d? (^) = -^dt (^) + d?- (^) 
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go respectively to exp [«(/)] and exp [—?(/)] at the d = limit. 
The recurrence relations for n = 2, m = 3 are solved by 

41 = 1 . ^sV = . 4?+ = , VA; > 2 , eSU^) = 1 



•-23- ~ '-' ' '-23- ~ -■- ' '-23- ~ (- ) '-23- ~ aq ' '-23- ~ 1 10' ^23-1^^ ~ 2_^ '-lg+ ^ 

fc=0 



Because C[0, 0, arccosf]+iS'[0, 0, arccosf] = 4= the linear combination ^23 ('^) = 72'^23+('^) 
leads to the Kepler/Coulomb eigenfunction in the d = limit. 

In fact, we obtain in this limit the right Kepler/ Coulomb triplet with energy E2 — ^^"^° 
and angular momenta 1 -limrf_>.o/2i = limrf_>.o/22 = — ^^, liiiirf_s.o -^23 = 0-: 



'2 — 9h2 



lim ij^2^' {u, v) = lim tj^^^' {u)C2i (^) = ( "^ J ^ ^ ^^'' ^"^ = ^n ('^' '^) 



\ 2 
mo; \ 2ma . 



limV^52^(M,t;) = limr752HM)^^J(t;) = ( — - ) re ^'' e ^"^ = V'i°li(^,V5) 
lim4?(u,t;) = \imri§ (yii'S iv) = \~J^) \^~^^) ^"^"^ ^ ^S?('"'V') 

6 Further comments 

We end this long survey by thinking about further possible extensions of these ideas and 
calculations to other classical integrable systems. 

1. The immediate temptation is to address the Kepler/Coulomb problem in M^. The A^- 
dimensional supersymmetric Kepler/ Coulomb problem has already been developed and 
fully solved in |12]. In our approach it would be easily doable in the N = 0/N = 
1 and N = N/N = N — 1 sectors; the wave functions in L^(]R+ x S^"^) would be 
organized in irreducible representations of SO(A^ + 1) -the bound states- or SO(A^, 1) - 
the scattering states-. In both cases, the wave functions are Kummer confluent functions 
times the spherical harmonics in L^(S^~^). Of course, the structure of the supersymmetric 

N 



Hilbert space is more complicated; there are N + 1 sub-spaces of dimension , , 

0, 1, ■ ■ ■ A^, so that: J2k=o [ u ] ~ ^^- This means that in the fc-subspaces such that 
2 < k < N — 2 only the procedure proposed in j32] would be effective. 

The supersymmetric modified Euler/Coulomb problem in M^ does not pose more diffi- 
culties than those encountered in this paper because all the N — 2 additional variables 
are cyclic. 

2. It will also be interesting to build the supersymmetric extension of the Kepler/Coulomb 
problem constrained to a sphere. This problem in the non supersymmetric framework 
was addressed independently by Higgs and Pronko in |58] and [59]. The supersymmetric 
generalization will require to deal with all the subtleties of spinors living in non-fiat 
manifolds. The metric, the A^-bein, the spin connection, and the like will enter the 
supercharges one way or another to make the system more intricate, see e.g. [Hj where 
supersymmetric systems in curvilinear coordinates are constructed. 

The Euler problem considered on a sphere, see [60j, is also a fairly well known integrable 
system with applications in celestial mechanics. It seems promising and interesting to 
work out the corresponding supersymmetric extension. 
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3. Another important integrable system is the Neumann problem [6T]: a particle forced to 
move on a sphere under the action of attractive elastic forces. This has been a source 
of inspiration for treatises on dynamical integrable systems, see |62j and [HS], and it has 
been applied, in the repulsive case, by some of us to study topological defects in non- 
linear sigma models with quadratic [64J and quartic [65] potentials. We believe that the 
supersymmetric extension of the quantum version of the Neumann problem will provide 
a physical example of the systems envisaged by Witten in [5] to construct a quantum 
derivation of Morse theory. 



',(0)/ 



(0), 



4. The Bosonic zero modes iIjqq {r, 0) in (3) and ipQ-{{u, v) in (19) have been easy to find. The 
Fermionic zero mode ipQ (a;i,a;2) (28) in the two center problem was discovered by means 
of supercharges defined in elliptic coordinates, translated back to Cartesian coordinates, 
and checked within a Mathematica environtment. A Fermionic zero mode in the one 
center problem, however, does not exist [2], [12]. This is very intriguing and compels us 
to study SUSY quantum mechanics in curvilinear coordinates in a more profound way. 
The way forward to and backward from the curvilinear to Cartesian coordinates of all 
these structures is highly non-trivial, see [H] for early attempts in this program. We think 
that the difficulties with the zero modes have a similar origin to the subtleties arising in 
the definition of spinors on curved manifolds. We plan to analyze this issue in a future 
publication. 
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